HEAT KERNEL ESTIMATES AND RELATED 
INEQUALITIES ON METRIC GRAPHS 
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Abstract. We consider metric graphs with KirchhofF boundary condi- 
tions. We study the intrinsic metric, volume doubling and a Poincare 
inequality. This enables us to prove a parabolic Harnack inequality. 
The proof involves various techniques from the theory of strongly lo- 
cal Dirichlet forms. Along our way we show Sobolev and Nash type 
inequalities and related heat kernel estimates. 



Introduction 

The aim of this article is to study solutions of the heat equation on a 
metric graph. A metric graph is by definition a combinatorial graph where 
the edges are considered as intervals, glued together according the combina- 
torial structure. Metric graphs have received lot of attention in recent years 
both of the point of view of mathematicians and application (see e.g. the 
conference proceedings [21 [8] and the survey [21 ) . 

Studying differential operators on metric graphs is similar to studying a 
system of differential expressions. However the combinatorial structure of 
the graph is closely related to the spectrum and the heat conductance on 
the metric graph, see for instance [Sj \T3\ [7]. Our approach is somewhat 
different: equipping the (minimal) Laplacian with the so-called standard 
boundary conditions allows us to treat the metric graph as a strongly local 
Dirichlet space. We first show the validity of a family of functional inequal- 
ities such as Sobolev and Nash type inequalities. The latter one is known to 
be equivalent to the notion of ultracontractivity of the associated semigroup 
generated by the Laplacian. This leads to the existence of the heat kernel 
and an upper Gaussian bound of the heat kernel (see [B]. This follows the 
work done in [3l |T6] and [TO] . 

The main result of this article will be a parabolic Harnack inequality, which 
is in fact equivalent to two sided Gaussian estimates of the heat kernel. This 
also improves the upper Gaussian bound mentioned above in the sense the 
the first one does not reflect the graph structure whereas the bound coming 
from the Harnack inequality does. A similar upper bound was already ob- 
tained in |20j in the case of compact graphs, i.e. graphs where the edge set is 
finite. In [11] the heat kernel was explicitly given via a combinatorial series 
expansion, but the class of graphs is restricted to generalized star graphs, 
i.e. compact graphs with a finite number of rays attached. Our result will 
generalize both in the sense that we allow a priori general infinite graphs to 
obtain asymptotic estimates for the heat kernel. Globally, as we will see, the 
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underlying combinatorial structure will govern the heat conductance, but lo- 
cally the graph behaves as in one dimension. Note that in [15] the question 
is adressed whether the parabolic Harnack inequality is stable for discrete 
graphs under changing the conductance matrices. Their method uses metric 
graph techniques, which is called cable system there. As a by-product they 
show that the parabolic Harnack inequality holds on the metric graph if and 
only if it holds on the discrete graph. Our approach is different here, since 
we emphasize on the metric graph while the combinatorial structure will not 
play a major role. Moreover the setting of Dirichlet spaces allows us to treat 
arbitrary positive weights on the edges, see subsection 11.31 
Our method goes back to the work of Moser ( [iTl [181 H] ) • These ideas were 
generalized to the setting of Riemannian manifolds by Saloff-Coste ( |24U23j ) 
and to the setting of strongly local Dirichlet spaces by Sturm (|27l[28]). We 
use these works to obtain the Harnack inequality. Part of our results goes 
back to the author's Diploma thesis from 2008. After this work was finished 
we learned about the recent article [22] and the preprint [1] which have some 
overlap with our results. 

The paper is organized as follows: in section 1 we introduce the notion of 
metric graphs and the standard Laplacian on it. In section 2 several analytic 
inequalities are proven, in particular a Sobolev inequality and a Poincare in- 
equality. The latter one is one half of the necessary input for validity of the 
Harnack inequality. The first one however helps us also to obtain a Nash 
inequality and associated heat kernel upper bounds in form of an ultracon- 
tractivity estimate of the heat semigroup. In section 3 we will have a further 
look at the local character of the geometry of the metric graph in terms of 
volume doubling and Poincare inequalities on balls. Hence this discussion 
leads then in section 4 to the statement of the Harnack inequality together 
with a rough sketch of the proof. Furthermore we state some of its impor- 
tant corollaries, among others we mention two sided Gaussian estimates for 
the heat kernel. Finally in the last section we have a look at some typical 
examples. 

1. Basic concepts and results 

1.1. Metric graphs. A metric graph consists, like combinatorial graphs, 
of a countable set V of vertices and of a countable set of edges E C 
V X V\{{v, v)\v G V}. In contrast to the discrete graph we want to consider 
the edges as intervals, glued together at the vertices. Let /:£'—)• (0, oo) 
the map equipping each edge with a (finite) length, i : E ^ V mapping 
each edge to its initial vertex and j : E ^ V mapping each edge to its 
terminal vertex according to the graph structure. We say that an edge E 
is incident to a vertex V if either i(e) = v or j{e) = v. In order to get 
a topological structure, we follow [3] and [10]: let X^. := {e} x (0, /(e)), 
Xy := V L) [J Xg and Xe := X^ U {i{e), j{e)}. Now consider the canonical 

homeomorphism X^ ^ (0, Z(e)), (e, t) t and extend it to a homeomor- 
phism TTe : Xe ^ [0,l{e)] such that 7re(i(e)) = and 7re(j(e)) = /(e). Note 
that this forces continuous functions u S C{Xr) to be continuous on each 
edge interior, i.e. the functions Ue ■= u o ir^^ are continuous and the limits 
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lim Ue{x) arc equal for all edges incident to v for all v & V. Furthermore 

we are able to define a metric structure on the graph. Fix x,y & and 
let p G be a polygon connecting x,y, i.e. po = x, pN = y and for all 
k G {1, .., — 1} there exists an edge such that pk,Pk+i ^ '^e^:■ Note that 
since we do not allow loops, this edge is unique and if we force the graph 
to be connected there is at least one polygon connecting x, y. The length of 
such a polygon is just the sum of all of its segments, i.e. 

N-l 

^(P) •= Yl kefc(Pfe+l) -7re,(pfe)| 

k=l 

and therefore we can define the metric as 

d{x,y) := inf{L(p) : p connects x and y}. 

To make sure that this defines a metric we assume from now on that the 
graph is connected and that the vertex degrees are finite for all vertices, that 
is 

dy := \{e e E : V e {i{e) , j {e)}}\ < oo for all v eV. 
In order to get a suitable measure on Xr, we define for Y C Xr 



J u{x) dm{x) := J u{x)dm{x). 



Y ^^^eDY 

where m is the the measure induced by the images of the Lebesgue measure 
on each (0, l{e)). The LP spaces are then given by 

LP{Xr,m) = ^LnO,l{e)). 

eeE 

For later purposes, we define the mean of a function u to be 

Uy '■= -j' u{x) dm{x) := — j u{x) dm{x). 

Y Y 

1.2. The energy form. Our next aim is to define the energy form. For 
this purpose we define for 1 < p < oo 

W'^iXr) := W''P{0, l{e)) n C{Xr) 

and 

W^'^iXr) ■.= W''P{Xr)nCo{Xr), 

where Co(^r) denotes the closure w.r.t. || • ||oo of the set of continuous 
functions with compact support. The energy form is then given by 

V = V{£) = W^^\Xr), 

lie) 



£{u,v) := XI / u'^{x)v'^{x) dx. 
eeEi 
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One can show (see for instance [14j) that this form is associated with the 
operator acting on functions as —-^ on each edge with domain 

LGCo(^r)n0Ty2,2(o^;(g)) . ^^/(^)^ol, 

I e&E e~D J 

where u'^{v) denotes the directional derivative towards v along e. This oper- 
ator is known as the Laplacian with Kirchhoff vertex conditions, which we 
will denote by A. 

Note that we implicitly impose Neumann boundary conditions on vertices 
with degree one. Furthermore although we explicitly do not allow loops, 
multiple edges and edges of infinite length, we implicitly do, since we can 
view any x G \ ^ as a vertex with degree 2. Functions in the domain of 
£ will then automatically fulfill the continuity condition. 
Note that, £^ is a regular and strongly local Dirichlet form with energy mea- 
sure 

dT{u{x)) = \u' {x)\'^ dm{x) = |ng(x)p dx 

(of. [3]). This allows us to define functions locally in the domain, i.e. 

u £ T>\oc -^^ u£ Lfoj, and Wip £ T> D Cc{S) : ipu £T>. 

Using the energy measure one can define the intrinsic metric g by 

Q{x,y) := sup{|n(x) - u{y)\ : u G W^^^{Xr) nC{Xr) and \u'\ < 1}. 

For a discussion of the intrinsic metric in the context of general strongly 
local Dirichlet forms see [26] . In our situation we have the following result: 

Proposition 1.1. For all x,y G X-p we have 

Q{x,y) = d{x,y). 

Proof. Let y £ Xr and dy{x) := d{x,y), resp. gy{x) := g{x,y). Clearly: 
dy £ C{Xy') and for x G Xg, we have 

dy{x) := mm{d{y,i(e)) + |7re(x) - 7re{i{e))\,d{y,j(e)) + |7re(2;) - 7re(i(e))|}, 

1 2 

which gives dT{dy) < dm and therefore dy G W^^'^ ('^r)- Now by definition 
of Q we have dy{x) < Qy{x) for all x,y £ Xr- 

Let now x,y £ X^ and u £ W^^^{Xr) n C{Xr) an arbitrary function with 
\u'\ < 1. Then we have 

\u{x)-u{y)\ = \UeiTTeix)) - UeiTreiy))\ 

< J K{t)\dt 

< |7re(x) - 7re(y)|, 

i.e. Q{x,y) < |7re(x) — 7re(?/)| on X^. Choose a path p with L(p) = d{x,y) + e 
for some e >0. This yields 

N-l N-l 

Q{x,y) < g{pk,pk+i) < l^efc(pfc) - vre^(pfc+i)| = d{x,y) +e, 

k=l k=l 
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where denotes the unique edge with pk,Pk+i C A'ej,. Now the claim 
follows, since e > was chosen arbitrary. □ 

1.3. More general models. We can also treat more general forms as fol- 
lows: let Ce : ^ R measurable for all e e E. Assume that there exists 
A > such that 

< A"^ < Ce{x) < A 
for all e G -E and x G X-p. Then we define the form 

V = V{£) = W^^\Xr), 

lie) 

£{u,v) :=''^^ / u'^{x)Vg{x) Ce{x)dx. 

Then, we can compare the energy measure dT of this form with the energy 
measure dT of £: for u G ^{E) = T>{£) we have 

A-^dr(u) < df{u) < AdT{u). 

This implies for instance the equivalence of the intrinsic metrics of both 
forms, i.e. for all x,y G Xr we have 

^~^Q{x,y) < d{x,y) < Ag{x,y). 

Moreover the Nash, the Poincare and the parabolic Harnack inequality, see 
sections 2.2, 3 and 4, together with their consequences hold true with con- 
stants depending additionally on A. We will show these results for the form 
£ and leave all modifications for the results on £ to the reader, as they are 
easily derived from the above mentioned comparison of the energy measures. 

2. One-dimensional inequalities 

2.1. Sobolev and Poincare type inequalities. The next few results re- 
flect the one dimensional behavior in terms of some useful Poincare and 
Sobolev-type inequalities. 

Note that the graph as metric space is compact if and only if the edge set 
is finite. 

Lemma 2.1. Let Xr be a compact metric graph and u G W^^'"'^(A'r). Then 
there exists c G Afp such that for all x G Xr and each G *P connecting c 
and X, we have 

\u{x)-u\< J \u{y)\dm{y) (X) 

where u = f u dm. 

Proof. Since u is continuous it attains its maximum and minimum, i.e. 

7i(.T-|_) = maxn and n(.x_) = minu. In particular — n < and 

— -u > 0. Without restriction let the inequalities be strict, otherwise 
we choose the point c below as x_ or x+ resp. Since the graph is connected, 
there is a path pxt connecting X- and x+. Note that u restricted to this 
path is continuous and changes the sign, hence there exists c G such 



that u(c) = u. Take now x G Xx^ arbitrary and a path p^. connectmg c and 
X. Then by the fundamental theorem of calculus we have 

\u{x) — u\ < j \u' {y)\ dm{y) . 



□ 

This result yields a whole family of Poincare and Sobolev-type inequali- 
ties. Recall that the diameter of a metric space is defined as 

diamA^r := sup{d(x,y) : x,y £ Xr}. 

Corollary 2.2. Let Xy he a compact metric graph and u G W^'^\Xy) for 
1 < p < oo . Then 

1—- - 

ll'U — < (diam^Vr) p | A'r] « ||u ||p 
for arbitrary 1 < p < oo and 1 < q < co. In particular 

\\u — u\\p < {diam. Xy) p\Xy\'p\\u\\p. 
Furthermore u fulfills the Sobolev inequality 

ll^lloo < I'YrI p + (diam ^Yr) ^ 11^ lip- 

Proof. Both families of inequalities follow directly from the previous lemma: 
By (U]) we get for all y e Xr 



\u{y) — u\ < J \u {y)\ dm{y) < \p\^ f \\u \\p < \ diam Xr\^ ^\\u\ 



Hence we have proven the first inequality for q = oo. For 1 < (7 < oo we 
then have 



\u{y) — u\'^ dm{y) j < \Xr\'' sup \u{y) — u\ 

y&Xr 



1I--1 I- 
< I diam Xy\ p | At | « 



1 \\u'\\p. 



This gives the desired inequality in the case 1 < g < 00. 

Now set q = 00 and use the (inverse) triangle inequality to obtain 



u 



00 ^ \^r\ ||u||i + I diam Afrl ^ \\u\\p 



which by the Holder inequality yields the claim. □ 

Corollary 2.3. Let X^ be a compact metric graph. Then the embedding 

W'^P{Xr) ^ CoiXr) 

is continuous for 1 < p < 00 and compact for 1 < p < 00. 

Proof. The continuity of the embedding is clear by the Sobolev inequality in 
the previous corollary. To show compactness, let (fn) be a bounded sequence 
in W^'P{Xr), i.e for all n E N we have ||/„||i,p < M for some M > 0. By 

||/«||oo < C||/n||l,p < cM 



the sequence is bounded for all x G Afp- Furthermore we have 

\fn{x) - fn{y)\ < I \f^{z)\ dm{z) < LipyJ (^j \f{z)f dm{z)j " 

Pi 

where p| is an arbitrary path connecting x and y and L(pi) denotes its 
length. Hence we get 

\fn{x) - Uy)\ < dix,yy^\\f\\i,p < Md{x,yy^, 

i.e. the sequence (/„) is equicontinuous. The Arzela-Ascoli theorem now 
gives the compactness. □ 

Note that the simplest compact metric graph is a closed interval. The 
results above are well known for closed intervals. So we have seen that 
compact metric graphs could be treated like closed intervals. The simplest 
non-compact metric graph is a half-line. The lemma below shows the anal- 
ogy of half-lines to non-compact metric graphs. 

Lemma 2.4. Let Xy he an infinite metric graph and u G W^''''{Xr) for 
1 < p < oo. Then we have 

\\u\\oc < {^)'^\\u\\l,p 

for 1 < p < oo, and 

\\u\\oo < Wu'Wl 

for p = 1. Let Y G Xr be open and pre-compact, then for all u G Wq'^{Y) 
we have 

and 

\W\\oo < \\u\\l- 

Proof. Fix X G Xp. Choose yn,Zn G Xr with d(yn,Zn) = const, and 
d{yn, x) — >■ oo. Then we have 

\uiyn) - uizn)\ < / \u'{z)\ dm{z). 



Vyn 



Since u G VF^'^(Ar) the right hand side tends to zero, and hence u{y) — > c 
for some constant c for d{x,y) oo. But since u G LJ'{Xr) it follows that 
c = 0. This implies 

\u{x)\ < ll^u'lli. 

For 1 < p < oo we apply this inequality to and by Holder inequality we 
get 



1. 1 ,„ ,P-i 



\u 



nl<P^\\uP-^u'\\f <p^\\u\\p^ <{P^) — {\\u\\p + \\u'\\p), 



where the last estimate follows from the well known inequality ab < + 
^b*^, for all a,b > 0. The second statement of the lemma follows form the 
first by extending u by zero on Xr. □ 
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Remark. As already mentioned the inequalities above reflect the one di- 
mensional structure of the graph. Once is established, our arguments 
proceed essentially as in the case of an interval or a semi-axis, see for instance 

m- 

2.2. Nash inequality and the heat kernel. The Sobolev estimates above 
are the starting point for us to deduce bounds on the heat kernel. The first 
step will be the following Nash type inequality. 

Lemma 2.5. Let be an infinite metric graph andu £ W^''^{Xr)riL^{Xr). 
Then there exists c > such that 

1 2 

\\u\\2 < clln'lll \\u\\i . 
Let Y C Xy pre- compact, then for all u G W^'^{Y) n L^{Y) we have 

1 ^ 

\\u\\2 < cdXrl^^Wny + ||n'||2)3 llullj*. 
Proof. Let x G Xp. Then Lemma 12.41 applied to gives 

|n(x)p < 2||W||i < 2||n||2||u'||2 
and hence ^ ^ 

< ^/2||m|||||u'||||m(x)| 
which yields by integrating over Xy 

1 1. 

||n||2 < \/2||u||| llu'llg 

so the first statement is proven. Analogous calculations starting with the 
Sobolev inequality 

||^i||oo < I'^rT'^IKIIi + Ik'lli- 
from Corollary 12.21 with p = 1 yield the second inequality. □ 

Using this inequality we can show ultracontractivity estimates of the as- 
sociated heat semigroup. This uses ideas going back to Nash [21j. We refer 
to [6] and references therein for further details and only sketch the proof. 

Theorem 2.6. There exists c > such that we have for all u £ L'^{Xr) and 
t > 

II —fA II -ill II 

||e n||oo < ct 4, \\u\\2- 
Moreover the heat semigroup e~'^* has an integral kernel p{t, x, y) for all 
t > which satisfies 

< p{t, X, y) < ct^2. 

Furthermore this kernel is C°^((0,oo) x X^. x Xf,i) for all e,e' G E, all even 
derivatives w.r.t. space are continuous and all odd derivatives w.r.t. space 
satisfy the Kirchhoff boundary conditions. 

Proof. The first statement follows from Theorem 2.4.6 [6j and the estimate 
on the heat kernel from Theorem 2.1.2. [6]. We only prove the last claim. 
We prove this along the lines of [|6j Theorem 5.2.1. 

By the spectral theorem for all / G L^(Afr) we have that e~*^/ G P(A"') 
for all n > 1 and t > 0. By Sobolev's inequality we have that for all e G .E 
and t > the mapping e~^^f\^ is C°°. Furthermore all even derivatives 
w.r.t. space are continuous and all odd derivatives w.r.t. space satisfy the 
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Kirchhoff b.c. in each vertex. Since the mapping t i— )• e~*^/ G D(A") is 
analytic, we can conclude that (t, x) ^ e-*'^/|^ is C°° for all e e E. For 
fixed t > and x E Afp the map / i— )• e~*'^/(x) is bounded on L^, hence 
there exists a{t, x) £ such that 

e-'^f{x) = {f,a{t,x)). 

This implies that {t, x) — t- a{t, x) is weakly infinitely differentiable on each 
edge interior and hence ([6]) also in norm. This implies norm continuity 
w.r.t. space for all even derivatives and the Kirchhoff b.c. for all odd ones. 
For g G we have by construction 

(e-*^/,5) = j {f,a{t,x))g{x) dm{x) 

and hence the identity 

e~*^g = j a{t,x)g{x) dm{x). 

Let now g,h £ then 

= J P{t,x,y)g{x)h{y) dm{x)dm{y) 

with the identity 

P{t,x,y) = (a(|,x),a(|,7/)) 

□ 

Now applying Theorem 2.2.3 [B] to the ultracontractivity estimate above, 
we arrive at a logarithmic Sobolev inequality, that is for all < n G 
W^'^iXr) n L^{Xr) we have 

I uHogfdm{x) < e£{u) + M{s)\\ff + log ||/||2 

for all e > and with M(e) = c — This implies Gaussian upper bounds 
for the heat kernel as shown in § 3.2 of [6], i.e. there exist constants a,b > 
such that we have 

1 ,g{x,y)^ 

< p{t, x, y) < at~ie~ « 

where q{x, y) is the intrinsic distance defined in section 1 and which coincides 
with the distance function d{x,y) by Proposition II. li 

Remark. In the case of compact graphs, the method of [H] was used in [2D] 
to obtain similar results. However the method there to establish a Nash 
inequality does not seem to work in the non-compact case. 
In explicit power series expansion for the heat kernel were obtained for 
graphs with finitely many vertices, i.e. for graphs where edges of infinite 
lengths are allowed. 
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3. Volume doubling and Poincare inequality 

In spite of one-dimensional behavior discussed in the previous section, the 
combinatorial structure causes some different local dimension effects in the 
sense of the following lemma. As we will see later, these result provide the 
necessary framework to use the results from [27] and [28]. 

Lemma 3.1 (Local volume doubling). Let x S and fix R > such that 
\B2Rix) n {v £ V : d^ > 2}\ < 1. Then for all < r < R we have 

< m{B2r{x)) < + 1) • m{Br{x)) < oo. 

In particular the volume doubling holds locally, that is, for all relative com- 
pact Y C X-p there exists v = viY) such that for all x G X-p, r > with 
B2r{x) CY we have 

< m{B2r{x)) < 2" ■ m{Br{x)) < oo. 

Proof. Without restriction \B2r{x) H {v G V : dy > 2}\ = 1, since the other 
case is trivial. Denote the vertex contained in this intersection by v. Then 
we have 

m{B2rix)) = m{Bd(x,v){x)) + m{B2r-d{v,x){v)) 
= 2d{v,x) + dy{2r - d{v,x)) 
= {2- dy)d{v,x) + 2rdy 

< (2 - dy)r + 2dyr 

< {l + ^)-m{Br{x)). 

The second claim follows easily from the first, since by relative compactness 
there exists a lower bound for the edge lengths and an upper bound for the 
vertex degrees. □ 

Theorem 3.2 (Local Poincare inequality). Let x £ X-p and fix R > 

such that \Bji{x) H {v e V : d^ > 2}\ < 1. Then for all < r < R and 
u G W^''^{Br{x)) we have 

\u{y) - Ub,(x)?' dm{y) < cpr"^ j \u {y)\^ dm{y) 

Br{x) Br{x) 

with cp = dy if V £ Bp{x) r\V or cp = 2 else. 

In particular the Poincare inequality holds locally, that is, for all relative 
compact Y C Xp there exists cp = cp(Y) such that for all x G Xp, r > 
with Br{x) C Y we have 

\uiy) - Ub,(x)\'^ dm{y) < cpr"^ J \u{y)\'^ dm{y) 

Br(x) Br{x) 

for alluG W^'^{Br{x)). 

Proof. Both inequalities follow from Corollary 12.21 with p = q = 2, since 
diamBr{x) = 2r and in the first case we have m{Br{x)) < dy ■ r. In general 
we have 

J \uiy) - UBr{x)\'^ dm{y) < m{Br{x))r J \u{y)\'^ dm{y). 

Br{x) Br{x) 
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But since Y is relatively compact the quantity cp := supj '"^^^^^^^ : B,.{x) C 
Y] is finite. □ 

If we assume a more uniform local structure, we can even prove stronger 
versions of the above statements. To do so, we need the following definition. 

Definition 3.3. Let Xy a metric graph as above. We say Xy is of hounded 
geometry if 

D := sup dy < oo 

and 

e := inf /(e) > 0. 

e&E 

As soon as the graph has a uniform lower bound on the edge lengths, we 
can also assume w.l.g. that the edge lengths are bounded uniformly from 
above by 2£ by introducing vertices of degree 2. 

For graphs of bounded geometry the volume doubling property holds uni- 
formly locally on Xr for all < r < | with uniform doubling constant 

CD = 1 + f . 

Furthermore for fixed L > | =: the doubling property then holds also 
uniformly for < r < L with a different constant depending on the choice 
of L. We briefly sketch the proof taken from |25j : it suffices to show that 
for all X £ X^ and all L > i' we have 

m{B^ix))<clm{BL{x)). (0) 

Consider therefore a maximal set B of points in B^_ii_{x) at distance at least 
4- apart, which forces the balls Bii_{y), y G 0, to be disjoint and contained 

2 

in Bii{x). Hence 

Y,m{B,^{y))<m{BL{x)). 

yeB 

Since B is maximal, the balls Bii{y) cover B ii_{x) and in particular the 

^ 2 

balls B2i'{y) cover B ei_{x). This yields 

m{B^^,_^{x)) < Y,m{B2fj{y)). 

Now by applying twice the doubling property we have ( [^ with a doubling 

constant of at most : choose n G N such that (n — 1)^ < L < n^. Then 
we have 

m{BL{x)) < clm{B._t^{x)) < ... 

^ 2 

< {clr-'miB^_^^_^^,{x)) < cfmiB,{x)). 
Summarizing the discussion above, we have proven the following corollary. 
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Corollary 3.4 (Uniform local volume doubling). Assume the graph is of 
bounded geometry and fix L > 0. Then there exists v = D, L) such that 
for all X G X-p and all r G (0, L) we have 

< m,{B2r{x)) < 2" ■ m{Br{x)) < oo. 

Remark. Associated with the doubling constant we define the local dimen- 
sion to be = ^"^1^2"^^ — 1- This is motivated by the fact that in R" we 
have vol{B2r{x)) = 2"'vol{Br{x)). So in this case the local dimension and 
the usual dimension agree. Note that every V > v defines a local dimension 
as well and that the local dimension may change from point to point. 
For further details concerning volume doubling and its consequences we refer 
to the book [25J. 

If the graph is of bounded geometry, then as consequence of Theorem [321 
the Poincare inequality holds true uniform locally as well: 

Corollary 3.5 (Uniform local Poincare inequality). Assume the graph is of 
hounded geometry and fix L > 0. Then there exists cp > such that for all 
< r < L and u G W'^''^{Br{x)) we have 

1 \u{y) - ub,{x)\^ dm{y) < cpr^ j \u {y)\ dm{y). 

Br(x) Br{x) 

Remark. Note that assuming the lower bound on the lengths forces the 
metric space (Afp, d), and hence (Afp, q), to be complete. 

4. The parabolic Harnack inequality 

The aim of this section is to state the Harnack inequality and to give 
a brief sketch of the ideas of its proof. These ideas go back to the funda- 
mental work of Moser [T7], [18] and [19] and were based on the validity of 
local Poincare and Sobolev inequalities. However in |23j it was shown that 
the doubling property and the Poincare inequality suffice to show the local 
Sobolev inequality. In our setting, i.e. in the case of strongly local Dirichlet 
forms this was generalized in |27j and |28j. So as we have shown volume 
doubling and the Poincare inequality in the previous sections, we obtain a 
Harnack inequality, see section 4.2 and for a brief sketch of its proof see 
section 4.3 as well. 

First we introduce the notion of weak solutions. 

4.1. Weak solutions. Let / C IR, be an open interval. We say that u : I — ?■ 
VF^'^(Afr), t^ut'is in L^(/ — > W^''^{Xy)) if u is Bochner measurable and 

j ll^t lll,2 < °0. 

Denote by W-^''^{Xy) the dual of W^''^{Xy) with respect to the inner product 
of L'^^Xy)., such that we have 

W^^'^{Xy) C L^{Xy) C W'^^'^{Xr). 

Since the embedding VF^'^(Afr) LP'{Xy) is dense and continuous, the same 
is true for the embedding of L'^[Xy) ^ W~'^''^[Xy). As usual we say that 
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V : I ^ W ^'^{Xr) is a generalized derivative of u : / — )• W^''^{Xr) if for all 

J u{t)Lp'{t) dt = - J v{t)Lp{t) dt 

I I 
and write f u for v. We set W{I x Xr) := {u G L'^{I Xr) : 3g G L'^{I 
W~^''^{Xr)) s.t. u' = g} with norm 

\\u\\w ■= {^j \\ut\\l,2 + 11^^*11-1,2 dt 
I 

We will deal with functions which are locally in W(/ x X^): Let G be an 
open subset of Xy and denote by Q = / x G the parabolic cylinder. Then 
y^\oc{Q) consists of all dm (8> dt- measurable functions on Q such that for 
all relatively compact G' C G and every relatively compact, open interval 
I' C I there exists a function u G W{I x X-p) with u = u on I' x G' . We 
say that a function u belongs to Wo(/ x G) if tt G VV(/ x Xp) and if for 
a.e. t £ I the function -u^ has compact support in G. Note that a function 
■u G yVo(/ X G) has only to vanish on / x dG, but neither on the upper 
nor on the lower boundary dl x G. Roughly speaking the space Wioc is the 
appropriate space of weak solutions and Wo the appropriate space of test 
functions with respect to the heat equation: 

Definition 4.1. We say that u is a local subsolution (resp. local supersolu- 
tion) of the equation 

—/\u = —u on Q 
dt 

if u € Wioc(Q) and 

£j{u, ^) := J £{u, ip)dt + j n, ^)dt<Q (*) 
J J 

(or £j{u,ip) > 0, resp.) for all relatively compact, open J C I and all 
nonnegative ip G Wo(Q). The function u is called a local solution if it is a 
local subsolution and a local supersolution. In this case ((*]) holds true with 
"=" for allipe Wo{Q). 

Remark. We have seen in section 3 that in the case of a graph of bounded 
geometry we have 

m{BL{x)) < Cj^ m{Bef_{x)). 

2 

for all -L > 0. Hence the volume growth is at most exponentially. This 
implies that 

oo 

dr = oo. 



1 



lnm(Br{x)) 



By Theorem 4 in [2Q\ this implies that £ is stochastically complete, which 
is equivalent to each one of the following: 

• Tfl = 1 for some (and then for all) t > 0, where Tt denotes the 
associated semigroup to £. 
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• For some (and then for all) a > every nonnegative solution u S 
L°°{Xt) of the equation (A + a)u = is identically 0. 

For more details on stochastical completeness see [Ml [S] . 

4.2. Harnack inequality. 

Theorem 4.2. Fix Q<£<r]<a<l and C, G (0, 1). Assume the metric 
graph is of hounded geometry and fix L > 0. Then there exists a constant 
ch such that for all balls Bj.{x), < r < L, x G and all s gM, we have 

sup u{s,y)<CH inf u{s,y) (PHI) 

whenever u is a nonnegative local solution of the equation —Au = in 
Q = (s,s + r^) X Bj.{x), where Q_ = (s + er^,s + r/r^) x (B(^r{x) and 
(5+ = (s + ar'^,s + r^) x B(^r{x)- Here the constant ch depends only on v, 
cp and additionally on the choice of parameters e, rj, a and but not on 
u, r, s and x. 

Of course, there are simple example of graphs which are not of bounded 
geometry, but then still a local version of the parabolic Harnack inequality 
is valid, see [2S]. 

Theorem 4.3. Fix < e < t] < a < 1, ( e (0,1) and let Y C Xr be 
precompact. Then there exists a constant ch = ch{Y) such that for all 
X E Xy, r > with Br{x) C Y and all s GM, we have 

sup u(t, y) < CH inf u(t, y) 
{t,y)&Q- {-t,y)&Q+ 

whenever u is a nonnegative local solution of the equation — Au = in 
Q = {s,s + r^) X Br{x), where Q- = (s + er^, s + ryr^) x B^rix) and = 
{s + ar'^ , s+r'^) X B(^rix)- Here the constant ch depends only onuiY), cpiY) 
and additionally on the choice of parameters e, tj, a and but not on u, r, 
s and X. 

4.3. The road to the parabolic Harnack inequality. In this section we 
give a survey on the proof of the Harnack inequality (|PHip in the case of a 
graph of bounded geometry as in theorem 14.21 

Assumption In what follows, we will assume that the graph is of 
bounded geometry, that is 

D := sup dy < CO 

and 

£ := inf /(e) > 0. 

e&E 

We have seen in section 3, that this implies that the volume doubling and 
the Poincare inequality hold locally uniform for all < r < L, where L > 
is fixed. 

As mentioned above, one crucial ingredient for the proof of (PHI) is a 
local version of the Sobolev inequality. It was proven by Saloff-Coste, see 
[23], that the volume doubling and Poincare inequality suffice to show it. 
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Theorem 4.4. There exists a constant cs > depending only on v and cp 
such that we have 

( / |„|S<i,„) ^ <,,.^( I |,.'(.)|^ / 

Br{x) Br{x) Br{x) 

for all < r < L, x £ Xy o,nd u G W^''^{Br{x)) with suppn C Br{x). 

Hence one has all tools to imitate the original proof by Moser to obtain the 
sub- and supersolution estimates. In the context of Riemannian manifolds 
this was done in [23]. In the case of strongly local Dirichlet spaces it was 
done in ^7] and yields in our setting: 

Theorem 4.5. • Fix 1 < p < oo and let 6 E (0, 1). Then, for any 

nonnegative local subsolution u of the equation —An = in Q = 
{s, s + r^) X Br{x) we have the estimate 



supuP <C6-''-'^{r'^m{Brix)))-^ // \u\Pdmdt, 

Qs Jj 

Q 

with Qs = {s + 6r'^, s + r^) x B(^i_sy{x), where the constant C > 
is independent of u, 6, s and the ball of radius < r < L. 

• Fix < p < oo and let 5 G (0,1). Then for any nonnegative local 
supersolution u of the equation —An = ^n inQ = (s, s + r^) x B^ix) 
we have the estimate 

s\XY>u-P <C5-''-^{r'^m{Br{x)))~^ If \u\-Pdmdt 
Qs JJ 

Q 

with Qs = {s + 6r'^, s + r^) x B^i_sy(x), where the constant C > 
is independent of u, 6, s and of the ball of radius < r < L. 

• Fix < po < 1 + f • Then for all < 6 < 1 and all < p < po 
any nonnegative local supersolution u of the equation —An = ^n in 
Q = {s,s + r^) X Br{x) satisfies 

\u\P<' drndty^' < (c5-^-'^ir'^m{Brix)))-^^ j \u\p dmdt, 

Q 

where Q'g = (s,s + (1 — (5)r^) x B(^i_syr{x)- Here, the constant C is 
independent of 5, p, u, s and of the ball of Radius < r < L but it 
depends on pQ . 

These results can be viewed as the easy part of the proof. The more 
technical one is to deduce from the supersolution estimates a weak Harnack 
inequality. A basic tool will be the following weighted Poincare inequality, 
which is once again a consequence of the volume doubling and the strong 
Poincare inequality, see |23] and j28] . 

Lemma 4.6. Fix 5 G (0, 1). Then there exists a constant cpi such that for 
all < r < L, xo e Xv and u £ W^''^{Brix)) 

J \u — ub^'^\'4^'^ dm < C'pr"^ j tp'^dT{u) 

Br{xo) Br(x) 
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where ii{y) = (1 - f^i^r^i^-iliifll)^ and ub,^ = J u^lP' dm/ J ■ip'^ dm. 

Br{x) Br{x) 

Here the constant C'^ depends only on Cp, v and 5. 

Using this inequality, one obtains estimates of the size of the logarithm 
of a local solution. We denote by m (8) A the product measure on Xy x R. 



u of the equation — Au = in Q = (s,s + r^) x Br{x), there exists a 
constant c = c{u, r) such that, for all n > 0, 



Lemma 4.7. Fix 6,t £ (0, 1). Then for any nonnegative local supersolution 

Br{x), there i 

{m(^ X){{{t,z) e K+ -Aogu < -fj,- c}) < Cr'^m{Br{x))fi^^ 

and 

{m(^ X){{{t,z) £ : logu > fi - c}) < Cr^m(5,.(x))^~^ 

where = {s + rr'^ , s + r'^) X B(^i_sy{x) and = {s, s + rr'^) x Bi^i_s)r{x). 
Here the constant C is independent of n > 0, u, s and of the ball of radius 
<r < L. 

Combining this estimate on logu with the supersolution estimates, one 
obtains via the so-called abstract lemma, see [25], the following weak Har- 
nack inequality. 

Theorem 4.8. Let 0<p<l + l,0<e<r]<a<l and < C < 1- Then 
any positive local supersolution u of —An = in Q = {s,s + r^) x Bj.{x), 
satisfies 

( Jl \u\P dmdt^ " < C{r'^m{Br{x)))p iniu, 
Q- 

where 

Q- = (s + er^, s + r/r^) x B{x, (r) and = (s + crr^, s + r^) x B{x, (r). 

Here the constant C is independent of u, s and of the ball Br{x) of radius 
<r < L. 

As mentioned already, this together with the subsolution estimate finally 
yields the Harnack inequality. 

4.4. Consequences. In this subsection we collect some well known conse- 
quences of the Harnack inequality to underline its power. Throughout we 
assume that the volume doubling and the Poincare inequality are satisfied 
local uniform. In the case that both properties are just satisfied locally, one 
has to change the assumptions in an obvious way. Note that in that case all 
intrinsic constants depend heavily on the subset Y. For the proofs we refer 
to chapter 5 of [25j. 

We start with the following Holder estimates for local solutions of the heat 
equation. 

Proposition 4.9. There exists constants a £ (0, 1) and C > such that 
for all < 2r < L and T £ H we have 



\u{s, y) — u{t, x)\ < C sup \u\ 



t\2 +d{y,z) 
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for all local solutions u of —An = on Q = [T — 4r^,T) x B2r{x) and 
{s,y),it,z) G {T-r^,T)xBr{x). 

If one considers time independent solutions one directly deduces the el- 
liptic Harnack inequality: 

Proposition 4.10. There exists a constant C > such that for all < 
2r < L we have 

sup u < C inf u 

Br(x) Br{x) 

for all nonnegative local solutions u of Au = on B2r{x). Furthermore if 
Y C A^r is relatively compact, then there exists a constant Cy > such that 

sup u < Cv inf u 

Y y 

for all nonnegative local solutions u of Au = onY. 

Proof. The first statement follows easily from the parabolic Harnack inequal- 
ity. The second statement follows by a chain of balls argument. □ 

Remark. For a different proof of the elliptic Harnack inequality see also 
[lOj . which uses the Moser iteration technique developed in [ITj. 

Clearly, by Proposition 14.91 also Holder estimates for solutions of the time 
independent equation hold true. 

Proposition 4.11. There exists constants a G (0,1) and C > such that 
for all < 2r < L we have 

\u{y) — u{z)\ < C sup \u\ 

B2r{x) 

for all nonnegative local solutions u of Au = on B2r{x) and y,z £ Br{x). 

Proposition 4.12. Assume now that the volume doubling and the Poincare 

inequality hold uniformly in r > on Xy- Then has the strong Liouville 

property, that is any solution of Au = on that is bounded from below (or 

above) is constant. Moreover, there exists an a > such that any solution 

u of Au = on At with the property 

1 , , 
lim — sup |n| = 

for some x G Xy, must be constant. 

The next consequence of Harnack's inequality improves the upper Gauss- 
ian bound of the heat kernel from section 2.2. Moreover, together with the 
stated lower bound, this estimates is in fact equivalent to the local uniform 
parabolic Harnack inequality and hence to the local uniform volume dou- 
bling and the Poincare inequality. For the proof we refer to ^27t i28j (see [23] 
as well). 

Proposition 4.13. There exists constants ci,C2,Ci,C2 > such that for 
all x,y £ Xy and all < t < we have 

Cl / d{x,yf\ /, N ^ C2 f d{x,yf^ 

■ exp TT— < P{t, X, y) < —7^ exp 



m(B^(x)) Cit 7 '^^-m(i?^(x)) C2t 
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Figure 1. A compact graph 
5. Examples 

In this section, we have a look at examples which appear frequently in 
the literature. 

5.1. Generalized star graphs. We call a graph a generalized star 
graph if we can decompose the edge set 

E = EfU\jEi 

where the set Ef is finite set of edges of finite length and the sets Ei give a 
finite set of the half-lines, see fig. 1 for a star center Ej. Hence the graph is of 
linear volume growth and the quantity sup "^(-^^(^)) appearing in Poincare's 

x,r 

inequality exists. Moreover, the Harnack inequality holds uniform in r > 
and X G Xr- These graphs where studied in [12], [13] and [TT] . 

5.2. The EucUdean lattice 'Z'^. As a second family of examples we con- 
sider the lattices Z'^, d G N. That is V = 'Z'^ and n ~ m if |n — m|i = 1, 
where | • |i denotes the £i-norm in IR*^. Moreover each edge has length one, 
see fig. 2. An easy calculation shows, that the volume doubling holds uni- 
form. Rather obvious is, that the Poincare inequality holds uniform as well. 
This was shown in [22j in a more general context of Cayley graphs associated 
with finitely generated groups of polynomial volume growth. 

Consider now the metric graph Z"', which has the same vertex and edge set. 
We then assign to every edge a length / S with > 0. We then 

define a quasi-isometry in a canonical way: denote by vrg the homeomorphism 
mapping Xg C Z"' to [0, 1], and denote by TTe the homeomorphism mapping 
Xe C to [0,/e]- Let Le : [0, 1] — ^ [0, /e], X ^ hx. The quasi-isometry J is 
then defined as Je = 7r~^oL~^o7fg. Hence y) <d{Jx,Jy) <l+d{x,y) 
and therefore we have Bi_j.{x) C J{Br{x)) C Bi^ri^). This implies that the 
volume doubling holds globally on Z'^, where the doubling constant depends 
now additionally on By the chain rule we see that the Poincare in- 

equality also holds uniform. Hence the Harnack inequality holds uniform 
and the asymptotic heat kernel behavior is the same as in Z'^. 

5.3. Trees. Typical examples of metric graphs which do not satisfy volume 
doubling or Poincare inequality globally are trees. We call a metric graph 
a tree, if it contains no cycles, i.e. there exists no continuous mapping 
7 : [0, 1] — )• Xr with 7(0) = 7(1) and which is injective on (0, 1) (see Fig. 3 for 
a finite tree) . Assume that the vertex degree is at least 2 for all vertices v £ V 
and that the edge lengths are equal to 1 . Then the volume is of exponentially 
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Figure 2. The metric graph 1? 




Figure 3. A tree 

growth. Moreover, the doubhng property cannot hold uniformly. If we 

assume further that the vertex degrees are uniformly bounded from above, 
then we know by section 3 that the doubling property and the Poincare 
inequality, and hence Harnack's inequality, hold locally uniform. 
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